A Common Universality Class for the Three— Dimensional Vortex Glass and Chiral 
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We present a Monte Carlo study of the d = 3 gauge glass and the XY-spin glass models in the 
vortex representation. We investigate the critical behavior of these models by a scaling analysis 
of the linear resistivity and current-voltage characteristics, both in the limits of zero and strong 
screening of the vortex-interactions. Without screening, both models show a glass transition at a 
finite temperature and, within the numerical accuracy, exhibit the same critical exponents: z ~ 3.1 
and v = 1.3±0.3. With strong screening, the finite temperature glass transition is destroyed in both 
cases and the same exponent v = 1.05 ± 0.1 is found at the resulting zero temperature transition. 



PACS numbers: 7550. Lk, 7540.Mg, 7460.-W, 0550. +q 



I. INTRODUCTION 

It has been sugg estedH that defects may collectively 
pin flux lines (vortices) in a type-II superconductor in 
a field, leading to a vortex glass phase with vanishing 
linear resistance. In many numerical studies of the vor- 
tex glass transition, a_ simple model called the "gauge 
glass" has been usedBlj A related model is the XY- 
spin glass, which has been studied extensively in order 
to understand the magnetic ordering of a variety of mag- 
netic compounds with random and frustrated interac- 
tions. XY-spin glasses are of special interest since they 
potentially exhibit two distinct kinds of ordering: spin 
glass ordering due to freezing of the spins, and "chiral 
glass" ordering due. to. freezing of local chiral (vottesd 
degrees of frccdomlrllS It has been well establishedll3H 
that the spins do not have a finite temperature spin g l a ss, 
transition in three dimensions, whereas KawamuraE-iO 
has argued that a finite temperature chiral glass transi- 
tion does occur. This intriguing claim provides one of 
the main motivations for the present study. 

In this paper we present a comprehensive Monte Carlo 
study of the vortex glass transition in the gauge glass 
model and the chiral glass transition in the XY-spin 
glass in three dimensions. We consider both the situ- 
ation where screening between the vortices is neglected, 
(which is the case in most of the earlier work) and also 
where there is strong screening of the vortices. We find 
that both with and without screening, the chiral glass 
and gauge glass have very similar behavior. Without 
screening they have a finite temperature transition with 
numerically very similar values for exponents, suggest- 
ing that they may lie in the same universality class. For 
both models, we find that screening destroys the finite 
temperature transition. 

Our paper is organized as follows: In Se c, |ll| we define 
the models under consideration. In Sec. UHpwe discuss 
the quantities that we calculate and explain the finite 



size scaling techniques used in the analysis. In Sec. IV 
we present our results for the gauge glass model without 
screening, and briefly review results for the gauge glass 



with screening that we found earlier in Ref. [ Q (referred 
to as WY). In Sec. we present results for the the XY- 
spin glass with and without screening. We summarize 
our results and draw our conclusions in Sec. 



VI 



II. THE MODELS 

In the absence of screening the Hamiltonian of both 
the XY-spin glass and the gauge glass can be written in 
the phase representation as 



H 



J cos( 



(1) 



where the <pi are interpreted either as phases of a su- 
perconducting order parameter (gauge glass) or as the 
angles of two-dimensional spins (spin glass). Here, J is 
the interaction strength (henceforth set to unity), and 
the sum is taken over all nearest neighbor sites on 
a simple cubic lattice. In the case of the gauge glass, the 
effects of the external magnetic field and the disorder are 
represented by quenched vector potentials Aij, taken to 
be uniformly distributed in the interval [0, 2ir]. In the 
case of the ±J XY-spin glass, the random sign of the 
bonds between spins is represented by quenched vector 
potentials taken randomly to be (+J) or ir (-J). 

The Hamiltonian (jl|) obviously possesses a U(l) sym- 
metry, i.e. the model is invariant under the transforma- 
tion cf>i — ► fa + C Vi, where C is a constant. For the 
gauge glass this is the only symmetry. However, for the 
XY-spin glass there is an additional "reflection" symme- 
try, 4>i -> -fc Mi. 

It is convenient to rewrite the Hamiltonian in such 
a way that the chiral (vortex) variables, which are our 
main concern, appear explicitly. This transformation 
involves replacing the cosine in Eq. ([!]) with the peri- 
odic Gaussian Villain function, separating spin wave and 
vortex variables ,j— and then performing fairly standard 
manipulationsE2rEJ to obtain 
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Here, the vortex variables n; S {0, ±1, ±2, ...} sit on the 
links of the dual lattice (which is also a simple cubic lat- 
tice here), G(i — j) is the lattice Green's function 



G(i-j) 



(2tt) S 



1 - cxp[i k ■ (r. t - rj)] 

L " &o 2En=l[l-COs(fe„)] ' 



E 



(3) 



(with d = 3), and the are quenched fluxes given 
by (1/2tt) times the directed sum of the quenched vec- 
tor potential Ay on the original lattice surrounding the 
link on the dual lattice on which lies. Due to peri- 
odic boundary conditions, we have the global constraints 
Si t>i = J2i n i = 0- There are also the local constraints, 
[V-n]t = [V-b]j = 0, where the latter just follows trivially 
from the definition of b^ as a lattice curl. 

Since the Hamiltonian only depends on — b^ it is 
convenient to discuss the distribution of the quenched, 
fluxes when all the weight is shifted into the intervalE3 
< bf < 1 (where a is a Cartesian component). For the 
gauge glass model, the distribution of the shifted b a is 
uniform, i.e. 



P(b a ) = 1 (0 < b a < 1) 
= (otherwise), 



(4) 



while for the ± J spin glass the shifted b a have a bi-modal 
distribution with equal weight at (corresponding to an 
unfrustrated square on the original lattice) and 1/2 (cor- 
responding to a frustrated square): 



P(b° 
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") + 5{b a 



(5) 



Recent work on the gauge glass model and on the XY- 
spin glass have investigated the role of screening of the 
vortex-vortex interactions, which is. .a relevant perturba- 
tion near the critical teraperaturaru. It was found by 
Wengel and Young (WY)0 and Ref. [ || , that the vortex 
glass phase vanishes when strong screening is included in 
the d = 3 gauge glass model, and subsequent work by 
Kawamura and Lo found the same effect for the chiral 
glass transition. We therefore also discuss the effects of 
screening here. 

In the vortex representation, the Hamiltonian is still 
represented by Eq. (|^) but now the interaction G(i — j) 
has the screened form 



G(i-j) 



{2n) 2 
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1 - exp[i k- (r» - r 3 )] 
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(6) 



where Ao is a bare screening length. Note that in the 
long wavelength limit, the denominator is just k 2 + Aq 2 . 

In the simulations presented here, we consider just two 
cases: (i) A = oo, where there is no screening and the in- 
teractions between the vortices are long range, and (ii) 
A — > 0, where there is strong screening. In the latter 
case G(r ^ 0) = (27rAo) 2 with corrections which are ex- 
ponentially small, i. e., of order exp(— r/Ao). Because 



J2i( n i ~~ k») = we can always add a constant to G(r) 
for all r without affecting the results. We therefore add 
— (27rAo) 2 , as a result of which the only interaction is on- 
site, and then divide the interaction by (27rAo) 2 to have a 
well defined limit for Ao — > 0. The resulting Hamiltonian 
then has the very simple form 



= i^( ni -b0 2 (Ao-0). 



(7) 



NoteEil, however, that TLy is not trivial because the local 
constraint [V • n]j = effectively generates interactions 
between the n^. 

To summarize, we study four models in this paper: 

1. The gauge glass in the absence of screening. The 
Hamiltonian is given by Eq. (|J) where the G(i — j) 
are given by 2a. (|3|), and the distribution of the 
fluxes (shiftedE3 into the interval from to 1) is 
given by Eq. (g). 

2. The gauge glass with strong screening. The Hamil- 
tonian is given by Eq. (^) in which the distribution 
of the (shifted) fluxes is given by Eq. (ji|) . 

3. The chiral glass (i.e. vortex degrees of freedom in 
the XY-spin glass) in the absence of screening. The 
Hamiltonian is given by Eq. (^J) where the G(i — j) 
are given by Eq. (|3|), and the distribution of the 
(shifted) fluxes is given by Eq. (||). 

4. The chiral glass with strong screening. The Hamil- 
tonian is given by Eq. (^) in which the distribution 
of the (shifted) fluxes is given by Eq. (JsJ) . 



III. DATA ANALYSIS 

We simulate the Hamiltonians in Eq. (||) and on 
simple cubic lattices with N = L 3 sites where 4 < L < 
12. Periodic boundary conditions are imposed. We start 
with configurations with all — 0, which clearly satisfies 
the constraints, and a Monte Carlo move consists of try- 
ing to create a loop of four vortices around a square. This 
trial state is accepted with probability l/(l+exp(/3A_E)), 
where AE is the change of energy and (3 = 1/T. Each 
time a loop is formed it generates a voltage AQ = ±1 
perpendicular to it's plane, the sign depending o# the 
orientation of the loop. This leads to a net voltage 



V(t) 



2e 



I V (t) 



with I v (t) 



LAt 



AQ(t), (8) 



where I v is the vortex-current and t denotes Monte 
Carlo "time" incremented by At for each attempted 
Monte Carlo move. We will work in units where h/(2e) — 
1, and we set At = 1/(3N) so that an attempt is made 
to create or destroy one vortex loop per square in each 
direction, on average, per unit time. 
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The linear resistivity can be calcula 
age fluctuations via the Kubo formula 



2d from the volt- 



Plh 



i £ At (V(t)V(0)). 



(9) 



Here, (• • •) denotes the combined thermal and disorder 
average. Near a second order phase transition the linear 
resistivity obeys the scaling law@ 



Ptt 



1 (T,L) = L-( 2 - d +^(L 1 /-( T _ Tc )) ; 



(10) 



where £ is the correlation length exponent, i.e., 

e~(r-r c )-", (n) 

z is the dynamical exponent, and p is a scaling function. 
At the critical temperature, p becomes a constant and 
therefore pi in (T c , L) ~ L-( 2 " d + z ). if we plot the ratio of 
Pii n for different system sizes against T, then 



ln[p lin (L)/p lin (L / )] 
ln[L/L'] 



= d-2- 



at T r 



(12) 



i.e., all curves for different pairs (L,L') should intersect 
and one can read off the values of T c and z. We will refer 
to this kind of data plot as the "intersection method" . 
With the values of T c and z determined by the intersec- 
tion method we can then use a scaling plot according to 
Eq. ( |l0| ) to obtain the value of v. 

In the case of strong screening we find a zero tem- 
perature transition and a plot according to Eq. (^) to 
determine v does not work, since z = oo because there 
is activated dynamical scaling at the T = transition. 
However, one can still obtain static exponents by mea- 
suring the voltage generated by a finite external current, 
i. e., by I—V characteristics. In real superconductors, 
transport currents generate a non-uniform magnetic field 
because of Ampere's law, V x B = J. It is inconvenient 
to simulate a non-uniform system, so instead we effec- 
tively assume that the current is the same everywhere so 
each vortex feels a Lorentz force x J. The scaling be- 
havior of the response to such a perturbation should be 
the same as that derived earlier for response to an actual 
transport currento. We can therefore use this approach 
to determine critical exponents, which is our objective. 
The Lorentz force biases the moves and sets up a net 
flow of vortices perpendicular to the current, whose time 
average gives the voltage according tcQ Eq. (|). 

To analyze our data we need to understand the scaling 
behavior of the J-V-curves near a-,second order phase 
transition. The scaling theory givesoO 



T- 



E 



d-1 



d-2 



V T 



(13) 



where E is the electric field, J the current density, r a 
relaxation time, and g is a scaling function. At a zero 
temperature transition one has 



S-T'", (14) 
so, in three dimensions, Eq. ( |l3|) becomes 

From this equation we can see that the current scale, 
Jnl, at which nonlinear behavior sets in varies with T as 
Jnl ~ T 1+2v . Since the linear resistivity is defined by 



Phn = lim -J, (16) 
and g(0) can be taken to be unity, we can write 



E 
Jp\h 



fJ 



J 



(17) 



Furthermore, we expect that near the T — transition, 
long time dynamics will be governed by activation over 
barriers. Hence we expect 

T 1+l > lin = - = Acxp(-A£(T)/T), (18) 
r 

where AE is the typical barrier that a vortex has to cross 
to move a distance £. One can define a barrier height 
exponent i/j by AE ~ ~ T~^ v in terms of which 



T 1+i > lin = Aexp(-C/T 1+ ^). 



(19) 



We are able to obtain a rough estimate for tp from our 
data of the linear resistivity. 

In a finite system, the I-V characteristics will also de- 
pend on the ratio L/t;. One can generalize the scaling 
function, Eq. (|l7|), to account for finite size effects as 
follows: 



E 
Jpih 



= 9 



J 



jn+2// 



L l/u T 



(20) 



Now we are left with a rather complicated scaling func- 
tion since it depends on two variables. To simplify the 
analysis we first estimate v by determining the current 
where E/(Jpn n ) = 2, at which point non-linear effects 
start to become significant. Denoting these values of J 
by Jnl, then, from Eq. (po|), it follows that 



Jnl 



= ~q[L 1 / u T 



(21) 



where g is another function. Hence we determine v by 
requiring that the scaling in Eq. (|2|) is satisfied. We 
then collect data for sizes and temperatures such that 
L 1 /"T is constant. The scaling function in Eq. ( p0| ) then 
only depends on one variable, and so data for E j J p\ m for 
different sizes should scale when plotted against J /T 1+2v , 
with the same value of v as obtained from the scaling of 
Jnl- We find, in fact, that the results are only weakly 
dependent on the second argument of Eq. (pp|) . 
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FIG. 1. Plot of ln[p lin (L)/p lin (L')]/ln[L/L'] versus T for 
the gauge glass with Ao — > oo. The curves intersect at 
T c — 0.93 ± 0.05. At the intersection point, the y-vahie is 
approximately —2.2, corresponding to zgg w 3.2. 

IV. RESULTS FOR THE GAUGE GLASS 

In this section we consider the critical behavior of the 
gauge glass model with and without screening. Recall 
that the distribution of the (shifted) fluxes is given by 
Eq. (@). 

A. No screening, Ao — > oo 

For the gauge glass with no screening we have mea- 
sured the linear resistivity py ln as a function of tempera- 
ture. In Fig. we show data of pu n plotted according to 
the intersection method vs. T for sizes L = 4,6,8. We 
were not able to include data from L = 10 into this plot 
since we could not equilibrate the systems down to the 
lowest temperatures (T = 0.8,0.9). All curves intersect 
at about T = 0.93 ± 0.05 indicating a phase transition 
to a vortex glass. The corresponding j/-axis value at the 
intersection point is 1 — z ks —2.2, therefore zgg ~ 3.2. 
Having established these values, we tried a scaling plot 
according to Eq. ([n]) and the result is shown in Fig. [|. 
Best scaling was achieved with T c = 0.93, zgg = 3 and 
j/qg = 1.3 ±0.3. Only far away from the transition point 
does one observe deviations from scaling, which is ex- 
pected for such small sizes and high temperatures, but 
the overall scaling works quite well. 

It is interesting to compare this result with earlier 
Monte Carlo simulations of the gauge glass without 
screening in the phase representation by Reger et al.Q 
These authors did a finite size scaling analysis of static 
quantities which indicated a finite temperature transi- 
tion, but they could not completely rule out the possibil- 
ity that the lower critical dimension is di ~ 3. The clear 
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L iA( T -T c ) 

FIG. 2. Scaling plot of the lineage resistivity for the gauge 
glass with Ao — > oo. Using T c — 0.93 and zgg = 3 from Fig. [I] 
we obtain vgg = 1-3 ± 0.3. 

intersection of the data in Fig. [I], however, strongly con- 
firms the notion that there is a finite temperature tran- 
sition in the three-dimensional gauge glass model, and 
hence di < 3. Additionally, our correlation length expo- 
neat z^gg agrees well with the estimate given by Reger et 
al.Q, ^gg = 1.3 ± 0.4. There is, however, a considerable 
difference between our estimate of the dynamic critical 
exponent z « 3.1 and theirs, z — 4.7 ± 0.7. It is pos- 
sible, though, that the dynamical universality classes of 
the models in the phase and vortex representations may 
be different, even though the static behavior is the same. 
If so, there is no contradiction in the results. 



B. Strong Screening, Ao —* 

In this paragraph we review quickly the results for the 
gauge glass model with strong screening found earlier 
by WY, in order to compare them in the next section 
with our data for the d = 3 XY-spin glass model with 
screening. As shown by WY the vortex glass transition 
in the gauge glass is destroyed by screening of the vortex- 
interactions. The main indication for the lack of a tran- 
sition at finite T was the absence of an intersection if 
the resistivity was plotted according to the intersection 
method. A scaling plot of the current -voltage character- 
istics for different temperatures and sizes also revealed 
T c = and uqq = 1.05 ± 0.1. Finally, the barrier ex- 
ponent ip, as defined in Eq. ([l^), was determined to be 
close to zero, so the conclusion was drawn that energy 
barriers diverge only weakly, possibly logarithmically, as 
one approaches the zero temperature transition. 
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FIG. 3. Plot of ln[p lin (L)/p lin (L')]/ln[L/L'] versus T for 
the ±J XY-spin glass and Ao — * oo. The curves intersect at 
T c = 0.94 ± 0.05. At the intersection point, the j/-value is 
approximately -2.2, corresponding to zcg — 3.2. 

V. RESULTS FOR THE ±J XY-SPIN GLASS 

Recall that the only difference between the ±J XY- 
spin glass and the gauge glass discussed in the last section 
is that the distribution of shifted fluxes is given by Eq. (^) 
rather than by Eq. (Q) . 

A. No screening, Ao — > oo 

As already discussed in the introduction, the ±J XY- 
spin glass is known to have no finite-temperatuEfi^tran- 
sition to an ordered state below four dimensionstll. For 
the d = 3, ± J model one can, however, identify a chiral 
glass transition in Monte Carlo simulations due to freez- 
ing out of the discrete degrees of freedom, as has been 
done by Kawamura et al.E-3 The associated chiral glass ex- 
ponents estimated in the piiase representation with pe- 
riodic boundary conditions^ are vcg — 1-5 ± 0.3 and 
t]cg = ~ 0r4 = ' = 0-2- Subsequent work with free boundary 
conditional^ finds similar values, isqg = 1-3 ± 0.2 and 
?7cg = -0.2 ±0.2. 

Figure || displays a plot of our data for pn n accord- 
ing to the intersection method vs. T for the XY-spin 
glass. One observes, very similarly to Fig. |l|, an inter- 
section point at T — 0.93 ± 0.05 and a dynamic critical 
exponent zcg ~ 3.2. Also, the scaling plot of pu n shows 
best results with almost the same values as in the long 
range gauge glass case, namely T c — 0.93, zqg = 3.1 
and i>cg = 1-3 ± 0.3. This result indicates a finite 
temperature transition into a chiral glass state for the 
d = 3 XY-spin glass and thereby confirms Monte Carlo 
results performed in the phase representations Very sur- 
prisingly, we find that our data for the linear resistivity is 
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FIG. 4. Scaling plot of the lineage resistivity for the ±J 
XY-spin glass with Ao — > oo. Using T c — 0.94 and zcg =3.1 
from Fig. Q we obtain vcg = 1.3 ± 0.3. 

virtually indistinguishable from the corresponding mea- 
surements of the gauge glass model. We observe a max- 
imum deviations of 1.5er. We will come back to this in 
the last section. 



B. Strong screening, Ao —* 

In Fig. H we show the linear resistivity plotted accord- 
ing to the intersection method vs. T. On can see, that 
there is no apparent intersection over the entire temper- 
ature range that we have been able to simulate, i. e., 
down to T = 0.07 for L < 8 and T = 0.1 for L < 12. 
At high temperatures all curves merge, since the correla- 
tion length becomes shorter then the system size and the 
data of pun for different sizes are the same. This rules 
out a transition down to 1/5 of the critical temperature 
of the system without disorder, T c = 0.331 (see WY), 
and therefore strongly suggests the absence of a chiral 
glass transition at finite terpaerature, in agreement with 
work by Kawamura and LiJl3 

Next we studied the current -voltage characteristics of 
our model in order to determine v. Figure ^| shows a scal- 
ing plot of different I-V curves according to Eq. J2C|) . 
From the scaling of the nonlinear current Jnl we esti- 
mated vcg — 1) and then chose sizes and temperatures 
for the data in Fig. || such that ^V^ccy = cons t.^ and 
hence the second argument in Eq. (|2^), remained roughly 
constant. The data is seen to scale very well with T c = 
and vcg — 1-05 ± 0.1. We also attempted scaling our 
data with an appropriate scaling function for finite T c , 
and found that scaling works only moderately well with 
T c = 0.04 and vcg = 1-05. We, therefore, conclude that 
the transition is very likely to occur at T c = 0, but we 
cannot completely rule out a finite, though extremely 
small T c . 
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FIG. 5. Plot of \n[p hn (L)/p lin (L')}/\n[L/L'] vs. T for the 
XY-spin glass model with Ao — > 0. In contrast to Fig. |^ 
there is no intersection over the entire temperature range, 
indicating the absence of a phase transition into a chiral glass. 



We also determined the barrier exponent ip by plotting 
T 2 pii n over 1/T as was done in Fig. 4 of WY for the gauge 
glass. The data for L = 12 follows almost a straight line 
indicating Arrhenius behavior and therefore ip ~ 0. As in 
the gauge glass case one has to be careful though, since 
such an estimate does not allow for finite-size corrections 
and is only observed over a small range of temperatures. 
It is also possible, that we only measure an effective ex- 
ponent and the true value of ip changes as one gets closer 
to T = 0. In any case, ip ~ would suggest, that barriers 
increase only very slowly, possibly logarithmically, as one 
approaches the zero temperature chiral glass transition. 

Again it is interesting to compare these results with 
those obtained by WY for the gauge glass with screening: 
they agree perfectly with in the errors, nam ely T c = 0, 
1.05 and ip 



FIG. 6. Scaling plot of the I—V characteristics with T c — 
and ucg = 1-05 ± 0.1, according to Eq. (poj), choosing sizes 
and temperatures such that L 1 ^ l ' CG T is roughly constant. 



vqq = l.ua and ip ~ 0, as described in Sec. [IV B . Not 
only do the final estimates of the exponents for the gauge 
glass and ± J XY-spin glass agree but also, as in the case 
without screening, the individual numerical values of pn n 
and data from the I—V characteristics all agree within 
the errorbars, the maximum discrepancy being 1.5a. 



VI. SUMMARY AND DISCUSSION 

In this article we have presented a Monte Carlo study 
of the gauge glass model and the chiral glass transition 
in the XY-spin glass model with and without screen- 
ing, in the vortex representation. We have computed 
dynamic quantities such as the linear resistivity and 
current-voltage characteristics and used finite-size scal- 
ing techniques to extract the critical behavior of these 
models. 

Our main results are: 



Screening 


T c 


V 


z 


i> 


Ao — ► oo 


0.94 ± 0.05 


1.3 ± 0.3 


« 3.1 


n/a 


Ao -> 





1.05 ± 0.1 


oo 


~ 



TABLE I. Critical temperatures and exponents of the pre- 
sumed common universality class of the gauge glass and the 
chiral glass in d — 3. v is the correlation length exponent, z 
is the dynamical exponent, and ip is the barrier exponent for 
the T — transition in the strong screening limit. 



1. In the absence of screening there is a finite tem- 
perature transition in both cases with numerically 
indistinguishable exponents given in Table. 1. 

2. In the presence of strong screening, there is a tran- 
sition at zero temperature in both cases. The cor- 
relation length exponent is the same for the two 
models, as shown in Table 1. 

3. Not only do the gauge glass transition and the chi- 
ral transition in the XY-spin glass model appear 
numerically to be in the same universality class, 
but even the individual data points for the current 
voltage characteristics are virtually indistinguish- 
able. 

Earlier work which provide d evi dence for a finite tem- 
perature chiral glass transitionli-alLa used the phase repre- 
sentation and constructed the chiralities (vortices) indi- 
rectly from the spin configurations. This is only sensible 
at moderate to low temperatures where correlations in 
the angles of nearest neighbor spins become significant. 
Our work is the first which demonstrates the existence of 
a chiral glass transition using the vortex representation. 
In our model vortices are well defined at all temperatures 



G 



and so we expect that the region over which scaling be- 
havior is obtained will be larger than in the earlier work 
in the phase representation. We therefore feel that our 
results make the existence of the chiral glass transition 
more convincing. 

There is also support in two dimensions for the idea 
that the chiral glass and gauge glass transitions are 
in the. [Same—universality class since in both cases one 
findiaOBfle T c = and v a 2. 

However it is unclear to us theoretically why the gauge 
glass and the chiral glass transition in the ±J XY-spin 
glass should be in the same universality class. For the 
XY-spin glass, the important low energy states are those 
where nf — bf = on links where bf = (corresponding 
to an unfrustrated square on the original lattice) and 
nf - bf = ±1/2 on links where bf = 1/2. Thus, as first 
noted by Villainja one has a random Ising model with 
long range antiferromagnetic interactions, 



H 



(22) 



where the ej are quenched variables taking values or 
1, and the Si are statistical Ising-like variables which 
take values ±1. For the gauge glass one cannot make 
an analogous transformation and Eq. (|^) corresponds to 
an Ising model in a random field, which is not expected 
to be in the same universality class as Eq. fl22]). Wc 
do not, therefore, understand why the numerical values 
of the exponents are the same within the uncertainties. 
Even more surprising is that the individual I-V values 
for the two models are virtually indistinguishable. We 
would expect there to be a more clearly visible difference 
in other properties. Perhaps for some reason, the random 
field aspect of the gauge glass is irrelevant, or perhaps 
the critical behaviors of the two models just happen by 
coincidence to be very close. It would be interesting to 
check our results by studying both models by alternative 
techniques such as domain wall renormalization group 
methods. 

The correlation length exponent for the unscreened 
models is also very similar to that of the Ising spin glassEEl 
with short range interactions. Again, it is not obvious 
to us why this should be the case. While the model in 
Eq. (|2|) has Ising variables, and the ingredients of ran- 
domness and frustration necessary for a spin glass, it also 
has long range interactions, unlike the Ising spin glass. 

Finally, it is noteworthy, that earlier results for the 
gauge glassO indicated that the universality class changes 
(and hence T c becomes zero) for any non-infinite value of 
tha.bare screening length. By contrast, Kawamura and 
LO have argued that the transition in the chiral glass 
persists down to a finite value of Ao- It would, therefore, 
also be interesting to study these models with an inter- 
mediate range of screening in the vortex representation. 
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